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(— i Abstract 

We argue that the existence of the cold dark matter is explained by primor- 
dial black holes. We show that a significant number of primordial black holes can 

q be formed in an axion-like curvaton model, in which the highly blue-tilted power 

spectrum of primordial curvature perturbations is achieved. It is found that the 
produced black holes with masses ~ 10 20 — 10 38 g account for the present cold dark 
matter. We also argue the possibility of forming the primordial black holes with 

q mass ~ 1O 5 M as seeds of the supermassive black holes. 
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1 Introduction 



The current cosmic microwave back ground (CMB) observations have revealed that the 
our present universe is filled with the unknown matter called dark matter, which cannot be 
explained within the framework of the well-established standard model of particle physics. 
The observed density parameter for the cold dark matter (CDM) is found by the WMAP 
[1] to be 

ttcDMh 2 = 0.1126 ±0.0036, (1) 

where h is the dimensionless Hubble parameter defined via the present Hubble parameter: 
H = 100/i km sec -1 Mpc -1 . In order to detect the dark matter, many experiments have 
been performed by now, but we have not found any meaningful signature yet. Therefore 
it is one of the most important problems of modern cosmology and particle physics to 
answer what the dark matter is. 

It is usually assumed that the dark matter is the weakly-interacting massive particles 
(WIMPs). The supersymmetric (SUSY) model [2], which is one of the most promising 
model beyond the standard model, naturally provides such WIMPs as the lightest su- 
persymmetric particle (LSP). Another promising candidate of dark matter is the axion, 
which is originally introduced to solve the strong CP problem in the standard model [3]. 
However, even if the SUSY particles or axion exist in the present universe, it may not 
be enough to explain the observed dark matter abundance depending on the model pa- 
rameters. In such a case, we are forced to demand another candidate for dark matter. 
It is known that the primordial black holes (PBHs), the black holes formed in the early 
universe [4], can behave like CDM. In this paper, therefore, we argue the scenario in which 
the currently observed abundance of CDM is explained by PBHs. 

PBHs are expected to be formed through the collapses of the high density regions 
caused by the large primordial density perturbations [5]. Light PBHs with mass smaller 
than 10 15 g are evaporated by now through the Hawking radiation [6], implying that 
only the PBHs with masses Mbh > 10 15 g can survive and contribute to the CDM. 
Furthermore, various cosmological and astrophysical constraints show that only the PBHs 
with mass 10 17 g < Mbh < 10 27 g can be the dominant component of the current CDM 
[7]. Although it is not easy to build the model in which a significant number of PBHs are 
formed, various models were proposed in the literature. Focusing on the inflation models, 
for examples, PBH formation was proposed in double inflation models [8, 9, 10, 11] or 
running mass inflation models [12, 13]. 

Another motivation to consider PBHs is the existence of supermassive black holes 
(SMBHs) at the center of galaxies [14, 15]. The observation of quasars (QSO) reveals 
that the SMBHs with mass M BH ~ 10 9 M exist at the redshift z fa 6 [16]. These black 
holes cannot be explained within the purely astrophysical mechanism, so we must rely 
on the primordial origin. If PBHs with sufficiently large mass, Mbh ^ 1O 3 M , can be 
formed in the early universe, they can play a role of the seed of SMBHs [17, 18]. 

In order for PBHs to form through the primordial density perturbations, we need 
the strongly blue-tilted power spectrum of the curvature perturbations, which gives the 
large density perturbations at small scales while the large scale density perturbations are 
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consistent with the CMB observation. However, the observation indicates that the scale 
dependence of the power spectrum is slightly red-tilted at large scales. This inconsistency 
is solved by employing a curvaton. The curvaton is originally introduced to generate the 
primordial large scale curvature perturbations instead of the inflaton [19]. In the curvaton 
model, a scalar field (called curvaton) acquires fluctuations during inflation and after 
inflation it decays into the standard model particles producing the adiabatic perturbation 
in the radiation dominated universe. 

In this paper, we consider that the curvaton is responsible for generating only the 
small-scale curvature perturbations while the large-scale perturbations are generated by 
the inflaton. After the decay of the curvaton, a significant number of PBHs can be 
formed through large density perturbations due to the curvaton. A specific model for the 
PBH formation with curvaton was proposed in [20], where three scalar fields (including 
inflaton and curvaton) with ad hoc couplings among them evolve non-trivially during 
inflation and leads to large density perturbations at small scales. Our mechanism for the 
PBH formation is completely different from that in Ref. [20]. We consider an axion-like 
curvaton field whose nature is very crucial for the PBH formation. Furthermore, axion- 
like fields often appear in various particle physics theories. We consider that one of such 
axion-like (curvaton) fields may play an important role for the PBH formation as studied 
in this paper. Ref. [12] also discussed PBH formation in curvaton model without concrete 
models. 

The remainder of the paper is organized as follows. In section 2, we show an axion-like 
curvaton model and we see the largely blue-tilted spectrum for the curvature perturbations 
is achieved. In section 3, we consider the PBH formation within the axion-like curvaton 
model. Section 4 is devoted to the conclusion. 

2 The axion-like curvaton model 
2.1 The potential of the curvaton 

In this section, we describe the axion-like curvaton model which was originally introduced 
in [21] (see also [22]) as an axion model with extremely blue-tilted spectrum of the isocur- 
vature perturbations. The model is built in the framework of supersymmetry and has the 
following superpotential: 

W = hS($$ - f), (2) 

where $ and S are chiral superfields whose i?-charges are +1, —1 and +2 respectively, 
and / is some energy scale and h is a dimensionless coupling constant. Here we assume 
that the model has a global U(l) symmetry and $, $ and S have charges +1, —1 and 0, 
respectively. In the limit of the glabal SUSY, the scalar potential is derived from (2) as 

V = h 2 \$$ - f\ 2 + h 2 \S\ 2 (\<$>\ 2 + |$| 2 ), (3) 

where, the scalar components are denoted by the same symbols as the superfields. Pro- 
vided that \S\ < / is satisfied, S tends to the origin and $ and 4> are settled on the flat 
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direction satisfying 

$$ = f with S = 0, (4) 

which makes the scalar potential (3) vanish. Hereafter, we assume that the flat condition 
(4) is always satisfied. Including the supergravity effects, the Hubble-induced mass terms 
are added to the scalar potential [23] as 

V H = Cl H 2 \<S>\ 2 + c 2 H 2 \$\ 2 + c s H 2 \S\ 2 , (5) 

where Ci, c 2 and c$ are numerical constants assumed to be real, positive and of order 
unity. In addition, there also exist the low energy SUSY breaking terms, 

V m = m 2 |$| 2 + m 2 |$| 2 + m||S| 2 , (6) 

where mi, m 2 and ms are soft masses of order of the gravitino mass. Here, because we 
are interested in the inflationary epoch, we neglect the low energy SUSY breaking mass 
terms. Thus, the flat direction is lifted by only the Hubble-induced mass terms (5) and 
the minimums of $ and $ are determined as 

f \ 1/4 / \ 1/4 

Now, we decompose the complex scalar fields into the radial and the phase components 

as 

$ = 7f^ exp ^+)' ^ = ^ ex p(^-)- ( 8 ) 

Then, along the flat direction, the massless direction is found as a linear combination of 
the phases, 9 = {9 + — 9_)/2. Without loss of generality, we can take <p ^> (p as the initial 
condition and neglect the dynamics of (p in the early epoch [22], so we follow the dynamics 
of only the complex scalar field $ = (pe l9 /a/2 whose potential is given by 

V v = \cHV. (9) 

Note that since the mass is comparable to the Hubble parameter, ip rolls down the po- 
tential somewhat rapidly during inflation. 

In our model, the curvaton is defined as the massless phase direction, which is expressed 
as a = ip6, so the curvaton field decreases with <p> until <p reaches the minimum y? m i n ~ / 
without changing 9. After that, the curvaton remains at the value a ~ f9 due to its 
masslessness. However, we assume that the U(l) symmetry is broken by some non- 
perturbative effect and a has the following potential in the low energy universe like the 
axion: 

'a 



V a = A 



cos 



\mW, (10) 



J 

where the second equality holds near the minimum <j m i n = and the curvaton mass is 
defined as m a = A 2 //. After the Hubble parameter becomes smaller than the curvaton 
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mass, the curvaton field starts to oscillate coherently with the initial amplitude cr, ~ f'6 
and behaves as matter. 

Let us derive the ratio r of the curvaton energy density to that of the radiation at 
the time of the curvaton decay. In order to estimate this, we consider the two cases : (i) 
the reheating occurs after the curvaton starts to oscillate and (ii) the curvaton starts to 
oscillate after the reheating. In case (i), which is denoted as m a > Y I (Yj : decay rate of 
the inflaton), we estimate r as 

P.faec) _ PM T r ^ 1/ f9\ 2 T R {Qim> Tl (n) 
' Pr{t dcc ) Pr(t R )T dcc 6\M P J T dcc ° rma ~M P ' 1 } 

and, on the other hand, in case (ii), m a < Tj, we obtain 

_ Pa(t osc ) _ 1 ( _f0_\ Tosc f < IK. do) 

~ Pr(t osc ) T dcc " 6 \M P ) T dcc M P ' [iZ) 

where the subscripts dec, osc and R correspond to the epochs of the curvaton decay, the 
curvaton oscillation and the reheating, respectively. Hereafter, we consider only the case 
of r < 1. 

The curvaton decays when the Hubble parameter becomes equal to the decay rate of 
the curvaton and the decay temperature of the curvaton is determined from the decay 
rate. Here we assume that the interaction of the curvaton with its decay product is 
suppressed by / like an axion, so we denote the decay rate of the curvaton as 

r -^f- < 13 > 

where k is a demensionless numerical constant assumed to be real, positive and smaller 
than 1. Then, the decay temperature of the curvaton is given by 

, x -1/4 

Tdec = - 5 (ifoJ ( r ^ 1/2 ' ( 14 ) 

where g* is the relativistic degrees of freedom. 



2.2 Generating the curvature perturbation 

To create a significant number of PBHs from the primordial density perturbations, we 
require the extremely blue spectrum with spectral index n s > 2 as we will see later. 
However, such a large spectral index is already ruled out by the CMB observation. In 
order not to contradict the observation, then, we build the model in which the almost 
scale-invariant large-scale curvature perturbations are generated by the inflaton and the 
small-scale curvature perturbations which are free from the CMB constraint are generated 
by the curvaton. Here, we investigate the possibility of PBH formation in the axion-like 
curvaton model introduced above. 
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The power spectrum of curvature perturbations is the sum of the contributions from 
the inflaton and the curvaton, which is written as 



V c (k) = P (M (k)+P <: , carv (k). (15) 

As mentioned above, the power spectrum is dominated by the first term in rhs. of (15) 
for small k and by the second term for large k. It is quite reasonable that the contribution 
from the inflaton is dominant until the perturbation scale at least k ~ 1 Mpc -1 leaves 
the horizon. For later convenience, we define k c as 

P(,cnrAkc) = 7W*c) - 2 X 1(T 9 (16) 

where we have used the CMB normalization [1] for the spectrum of the large-scale curva- 
ture perturbations. Using this definition, our requirement is denoted as 

P<,curv(*0 < V (M {k) ~ 2 x 1(T 9 for k < k c and k c > 1 Mpc" 1 . (17) 

The curvature perturbation spectrum from the curvaton is given by [24] 



Pc,curv(*0 [ A + 3r ) ( 2 n) a(ky (4 + 3r)(27rp(Jfe)0) (18) 

where the argument k entering in a and <p indicates the value when the scale k leaves the 
horizon. Note that, after (p reaches the minimum, the power spectrum takes the constant 
value given by 

*W(*) = « (jf^:) {^j§) for k>K > ^ 

where k* is defined as the scale leaving the horizon at the time ip reaches the minimum 
~ /. With use of the curvaton spectral index n a , the scale dependence of the power 
spectrum of the curvature perturbation is expressed as 

/ & \ n "- 1 

Pc,cu r v(k) = P Ccmv (k c )i-\ for k<K. (20) 

Combining this and the <p dependence, 

P C ,cnrv(k) = P C , cmv (k c ) (^) > (21) 



we obtain the relation 
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The spectral index of the curvaton is calculated by solving the equation of motion of 
ip with potential (9), 

<p + 3Hp + cH 2 ip = 0, (23) 
whose solution during inflation (H ~ const.) is given by 



tp oc e' xm oc k~ x with A=|-|i/l-^c. (24) 

2 2V 9 



Together with (20) and (21) the spectral index is given by 



n CT -l = 3-3^1 -^c, (25) 

so we can obtain the extremely blue spectrum such as n a ~ 2 - 4 with appropriate choice 
of c [21]. 

3 The PBH formation 

In this section, we consider the formation of PBHs in our model. It is well-known that 
PBHs can be formed by collapse of overdense regions in the radiation-dominated universe 
and their mass is as large as the horizon mass at the formation time [25, 26, 27], which is 
given by 



47T 



M BH = ^rP r H~ 3 ~ 0.05M G 



3 rr " ® ^100 y VGevy 

°Vioo; VMpc"V 



(26) 

-V6/ U. \- 2 



where p r represents the energy density of the radiation and the subscript / represents 
the time of the PBH formation. Here, we assume that PBHs cannot be formed from the 
curvaton itself, so the PBH formation takes place after the decay of the curvaton into 
radiation. PBHs with Mbh > 10 15 g do not evaporate through the Hawking radiation [6] 
until now and their abundance can contribute to the present CDM density. The current 
density parameter for such PBHs is calculated as 



PPBH,eq n w 1 n g R ( M ® 



n^^^^xiif;!-^)^, (27) 

where the subscript eq corresponds the time of the matter-radiation equality and (3 is 
defined as the energy density fraction of the PBHs at the PBH formation, which is denoted 
as = PYBnit f) / Ptotif f) ■ Various cosmological and astrophysical constraints are imposed 
on (3 [7], from which only the PBHs with mass ~ 10 17 g - 10 27 g can be the dominant 
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component of the dark matter. For such PBHs, the constraint on (3 comes from the 
current observational value of the dark matter density (1), which implies 



Assuming that the PBHs are created by the collapse of overdense regions with pri- 
mordial gaussian density perturbations, (3 is estimated as [28] 

^^) ex p(-18^))' ^ 

where 5 = Sp/p is the density contrast and (5 2 ) is its variance. The variance of the density 
perturbations is related to the power spectrum of curvature perturbations, and, in the 
comoving gauge in which the curvature perturbation is expressed as TZ, which coincides 
with C well outside the horizon, the following relation is known: 

4(1 _i_ w )2 

"•<*> = h^ mkl (30) 

at the time the scale k leaves the horizon [24] . w is determined from the relation between 
the pressure and the energy density of the cosmic fluid P = wp and it takes 1/3 in 
radiation dominated universe. The variance of the density perturbation smoothed over 
the scale R is estimated as 

(5\R)) = J W 2 (kR)V s (k)^, (31) 

where W(kR) is the window function in Fourier space. Assuming the Gaussian window 
function W(kR) = exp(— k 2 R 2 /2) and taking into account 7\ Viz and (19) and (20), 
we can approximate the variance as 

(S\R)) = ^W(k) [(hR)-^S((n c - l)/2, klR 2 ) + E^klR 2 )} , (32) 
where 7(0,2;) and E x are defined as 

j(a,x)= / t^e^dt, and E x {x) = / —dt. (33) 

If we wrire 

(6 2 (k- 1 ))=aV (: , CUIV (k) J (34) 

the numerical coefficient a is taken to be 0.1-4 as shown in Fig. 1. 

We show the energy density fraction of PBH in terms of V^ iCUVV in Fig. 2(a) in the case 
of a = 1 (solid red line) and a = 0.1 (dashed green line). The dotted blue line (the small- 
dotted magenta line) corresponds to the upper limit in the case of Mbh = 10 27 (10 17 ) g, 
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Figure 1: The ratio of the smoothed variance of density perturbation to the power 
spectrum of curvature perturbation are shown. The horizontal axis is the wave number 
corresponding to the smoothing scale as k = R~ l divided by fc*. This curve is independent 
of n„. 



which comes from the current observation of the CDM density. In order for PBHs to be 
the dominant component of dark matter, the required value of curvature perturbation is 
V Ccuiv ~2x KT 3 (2 x KT 2 ) for a = 1 (0.1). Substituting (18) and (34) into (29) and 
taking <p(k) ~ /, the constraint (28) is rewritten in terms of H- m f/(f6) shown in Fig. 2(b). 
In this figure, the thick (thin) solid red line corresponds to r = 1 and a = 1 (0.1) and the 
thick dashed green line corresponds to r = 0.1 and a — 1. The breaking point of each line 
corresponds to the point at whch the quantum fluctuation of the curvaton, 5a = H ini /2ir, 
becomes /. If H in {/2ir > f, the amplitude of the quantum fluctuations of S overtakes the 
critical value /, which invalidates our underlying assumption (4). From Fig. 2(b) we need 
r ~ 1 and f6 ~ ifj n f to account for the present dark matter abundance. 

Now let us estimate the mass spectrum of the PBHs in the present model. With 
R = k^ 1 and Eq.(26) we rewrite the smoothed variance (32) in terms of PBH masses as 



(5 2 (M BH )> = ^V Ccuvv (K 



Mi 



1311 



(n CT -l)/2 
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- 1 Mi 



BH 



M* 



M 



BH 



M* 



(35) 



where M* is the mass of PBH formed when the scale k* enters the horizon. Using (35), we 
can calculate the mass function, which is defined as the number of PBHs per comoving 
volume whose mass range is Mbh ~ Mbh + <^Mbh, as [29] : 



dn 



PBH 



dM, 



BH 



1 pr(Q ( M, 



18tt M 2 h 



M 



BH 



1/2 



d\n(5 2 (M BR )) 



din M BH 



/3(M 



BH J 



(<5 2 (M BH )> 



(36) 



where p r {t*) is the radiation energy per comoving volume at the time when the scale k* 
enters the horizon and /3(Mbh) is the density fraction of PBH whose mass is Mbh- Since 
we assume that PBHs are never formed before the curvaton decays, the mass spectrum 
has the lower cutoff M m i n , which corresponds to the mass of PBH formed just after the 
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Figure 2: The energy density fraction of the PBH at the formation is shown. The 
horizontal axis correspond to Pf jCurv in Fig. 2(a) and H ini / f9 in Fig. 2(b). In Fig. 2(a), 
the solid red line and the dashed green line correspond to a = 1 and a = 0.1 respectively. 
In Fig. 2(b), the thick (thin) solid red line corresponds to r = 1 and a = 1 (0.1) and the 
thick dashed green line corresponds to r = 0.1 and a — 1. Breaking point of each line in 
Fig. 2(b) corresponds to 5a /a = 1. The dotted blue line (the small-dotted magenta line) 
corresponds to the upper limit in the case of Mbh = 10 27 (10 17 ) g, which comes from the 
current observational value of the CDM density parameter : I^cdm = 0.23. 

curvaton decays. 1 The mass spectrum of PBH is shown in Fig. 3. The solid red and 
dashed green lines correspond to M min /M* = 10~ 8 and M m - m /M* = 10~ 3 respectively and 
they are normalized by the their own peak values. The mass spectrum depends only on 
M min and it is independent of n a . It is clear that the dominant contribution to the energy 
density of PBHs comes from the smaller mass PBHs, so the constraint on the initial PBH 
abundance should be applied to the PBHs with M min . In particular, for M-^^/M^ = 10~ 8 , 
it is seen that the number of PBHs with mass larger than ~ lO^M* decreases drastically. 
This is due to the sudden decreasing of a, which implies the sudden decreasing of (S 2 ), 
for k < lO 2 /^ (see Fig. 1). Thus we can obtain a very narrow mass spectrum by tuning 
M min /M* as 10~ 3 -10- 2 , which is required to explain SMBHs [30]. 

Now, we investigate the parameters allowing the formation of PBHs which eventually 
becomes the dominant component of the CDM. We also impose several conditions to build 
the viable scenario, which are listed below. 

• Going back to the time when the pivot scale k p = 0.002 Mpc^ 1 leaves the horizon, 
<p should be smaller than the Planck scale. Using (22) and taking into account 

1 Here we assume r < 1 which means the curvaton never dominate the universe before decaying and 
PBHs are never formed from the curvaton itself. If we allow r > 1 and the PBH formation from the 
curvaton (scalar coherent oscillation) itself, the minimum mass of PBH may be given by the horizon mass 
when the curvaton starts to dominate the universe. 
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M BH / M* 



Figure 3: The mass spectrum of PBH, dnpBK/ 'dM-en, is shown. The solid red line and 
dashed green line corresponds to M min /M* = 1CT 8 and M min /M* = 1CT 3 respectively and 
they are normalized by the their own peak values. They are independent of n a . 



f « H ini /6 and (p(k c ) « 10 3 /, we get 

/ L \ -K-l)/2 

H inf < 2 (^- 1 )/2io- 3 (^+ 1 )/^M P f^j . (37) 

Combining this with the constraint from the tensor-to-scalar ratio [1], ifj n f < 5 x 
l(T 5 Mp, we get 

H in{ < min [ 2 (nCT - 1)/2 10~ 3(nCT+1)/2 ^M P , 5 x 10~ 5 M P ] , (38) 
where we set k c = 1 Mpc -1 . 

The mass of PBHs formed when the scale k* reenter the horizon, M*, is calculated 
with use of (16), (20) and (26). Since M* is larger than the minimum mass of PBHs, 
we obtain the following condition: 

1/6 / , \ -2 /T) sx -2/(n CT -l) 



M mm < 1VU - I x 10 g^ oo j ^ Mpc _ x J ^ 2 x io _ 3 J 



(39) 



The curvaton should decay before the Big Bang nucleosynthesis (BBN), that is 
Tdec > 1 MeV. From (26), the minimum mass of PBHs is related to T^ ec as 

t — 1x103 Gev (i,)" 1/4 1/2 (4o) 
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Hence the minimum mass of PBHs is constrained as 




1/2 



M min < 1 x 10 38 g 



(41) 



which is combined with (39) and we get 




[ 2 x 10 



46-12/(^-1) 



g, lxl0 38 g], 



(42) 



where we set g* « 100, k c = 1 Mpc and V^ tCurv (K) = 2 x 10 3 . 

• The reheating temperature is also constrained. In the case of Tj < m a , from (11) 
and (40) we obtain 



where we take f6 ps H mi , r ~ 1 and g* ~ 100. In the case of Tj > m a , on 
the other hand, we get the similar relation by simply replacing Tr with T osc . The 
reheating temperature or the curvaton oscillation temperature is constrained from 
the inequality (38). 

We show the parameter space allowing for the PBH to take a role of the dominant 
component of the CDM in Fig. 4 in the case of Fj < m a and k c — 1 Mpc -1 . The allowed 
region is inside the respective contours. The dashed-and-dotted-cyan line is the lower 
limit on the PBH mass, which comes from the current upper limit on the tensor-to-scalar 
ratio. For the PBH to take a role of the dominant component of the CDM, we need the 
somewhat high reheating temperature T R > 10 12 GeV. 

Then, we investigate the allowed region of our model parameters, / and m a , in which 
the current dark matter density can be explained by the PBH. The allowed parameter 
region becomes much narrower than Fig. 4 if we take into account the decay rate formula 
(13). The relation between the decay temperature of the curvaton and PBH mass given 
by (40) is rewritten as 



which constrains the allowed region in / - m a plane when 10 g < M m - m < 10 g is 
imposed. In the case of m a > Tj, (11) is translated into the following inequality : 



where we set 6 = 1 and g* ~ 100. In the case of m a < Tj, on the other hand, the relation 
from (12) becomes approximately same as (45) but replacing > with ~. Moreover, the 




(43) 




(44) 




(45) 
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T R [GeV] 



Figure 4: The allowed parameter region in which PBHs become the dominant component 
of the CDM in our model is shown. The allowed region is inside the respective contours. 
The dotted-blue line and the solid-red line correspond to the boundary in the case of n a = 
3 and 2 respectively. The dashed-and-dotted-cyan line is the lower limit on the PBH mass 
coming from the upper limit on the tensor-to-scalar ratio. We have taken r = 1 and 6 = 1. 

constraints (38) is trivially translated into upper bound on / by simply replacing H in f 
with f9. The contraint (42) is rewritten as 

where M c is defined as the right-hand-side of (42). 

We summarize the above constraints in Fig. 5(a) and Fig. 5(b) for n a = 2. In the 
case of m a > T/, the arrowed regions are inside the thick solid-red (dashed-green) line 
for k = 1 (0.01) and, in the opposite case, m a < Tj, they are on the lower boundary of 
these lines. In addition, it must be below the thick (thin) dashed-and-dotted-cyan lines 
corresponding to the upper bound of the PBH mass, 10 g for k = 1 (0.01). From these, 
it is found that / and m a must be / ~ 5 x 10 13 - 10 14 GeV, ~ 5 x 10 5 - 10 8 GeV and 
A ~ 10 10 - 10 11 GeV to explain the current CDM abundance. 

Another outcome of our model is the possibility of explaining the seeds of SMBHs. 
The initial mass fraction of PBHs as seeds of SMBHs is constrained by the observed 
comoving number density of QSOs: a 3 riQso — (6 ± 2) x 10~ 10 Mpc -3 [16]. The comoving 
number density of PBHs is given by 

/ \ -1/4 / M \ 3/2 

aVBH^xiO'VMpc-^) (J£) , (47) 
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Figure 5: The allowed regions for the PBH to be the dominant dark matter in / - m a 
plane and / - A plane are shown. The allowed regions are inside the thick solid-red 
(dashed-green) lines and below the thick (thin) dashed-and-dotted-cyan lines for k = 
1 (0.01) in the case of m a > Tj. The dashed-and-dotted-cyan lines correspond to the 
upper bound of thePBH mass 10 27 g for k — 1 (0.01). We have taken n a = 2 and = 1 
in both figures. 



so, compared with a 3 nq SO , (3 is estimated as 

~ 2 x 10- 21 ( V4 ( -^-) V2 . (48) 

' V 10 / V 1O5M 0/ 

Since a quite large mass and a narrow mass spectrum of the PBH are needed to explain the 
SMBH, we set M min ~ M* which leads a ~ 0.1, H ml /f6 ~ 2 and V c , cmv (K) ~lx 10~ 2 . 
Then the parameter space is constrained by the same way as those of the PBH dark 
matter case and we summarize it in Fig. 6. 

In Fig. 6(a), the allowed region is inside the solid-red line (the dotted-blue line) for 
n a = 2.5 (n a = 3). The small-dotted-magenta line (the dashed-and-dotted-cyan line) 
corresponds to M BH = 1O 5 M (1O 4 M ), on which the SMBH is explained by the PBH. In 
Fig. 6(b), the allowed region is inside the solid-red line (the dashed-green line) for k = 1 
(0.01). The thick (thin) dashed-and-dotted-cyan line corresponds to M B h = 1O 5 M for 
k = 1 (0.01). We found that our model can provide the seeds of SMBHs for T R > 10 9 GeV, 
/ ~ 10 12 GeV, m a ~ 0.5 - 100 GeV and A ~ 10 6 - 10 7 GeV. 

However, the SMBHs cannot be a significant part of the dark matter density of the 
universe. Fortunately, various axion-like particles often appears in particle physics theo- 
ries. One of them may be the curvaton which is responsible for SMBHs as discussed above. 
Another axion field can play a role of the usual QCD axion which solves the strong CP 
problem. If the Peccei-Quinn scale f a is ~ 10 12 GeV, the QCD axion can account for the 
dark matter of the universe. The coincidence of two independent scales / ~ f a ~ 10 12 GeV 
may be very interesting. Furthermore, it is pointed out that the axion dark matter is a 
good candidate consistent with the presence of the primordial SMBHs [31]. The required 
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(a) (b) 

Figure 6: The allowed regions for the PBH to be the seed of the SMBHs in Mbh _ Tr plane 
and / - Trig plane are shown in Fig. 6(a) and Fig. 6(b) respectively. In Fig. 6(a), the allowed 
region is inside the solid-red line (the dotted-blue line) for n a = 2.5 (n a = 3) and the 
small-dotted-magenta line (dashed-and-dotted-cyan line) corresponds to Mbh = 1O 5 M 
(1O 4 M ). In Fig. 6(b), the allowed region is inside the solid-red line (the dashed-green) 
and the thick (thin) dashed-and-dotted-cyan line corresponds to Mbh = 1O 5 M for k — 1 
(0.01). We have taken r = 1 and 9 = 1 in both figures and n a = 3 in Fig. 6(b). 



scale A ~ 10 6 - 10 GeV is coincide with the SUSY breaking scale when it is mediated by 
gauge interactions, which suggests that the dynamics generates the curvaton mass may 
be related to physics of SUSY breaking. 



4 Conclusion 

We considered the axion-like curvaton model based on the SUSY, in which the curvaton is 
identified as the massless phase direction contained in some complex scalar field. Because 
of the Hubble-induced mass of the radial part of the complex scalar field, the power 
spectrum of the curvature perturbations from the curvaton becomes extremely blue such 
as n a = 2 - 4. In order not to contradict with the WMAP observation of the spectral 
index, in our model, the large scale perturbations (k < 1 Mpc -1 ) are generated by the 
inflaton giving the almost scale-invariant power spectrum and the contribution from the 
curvaton becomes significant at sufficiently small scales. We showed that, by use of such a 
extremely blue spectrum, the PBHs are formed from the collapse of the overdense regions 
and the produced PBHs have a rather broad mass spectrum. It is found that, in a certain 
parameter region, PBHs with mass 10 17 — 10 27 g can eventually become the dominant 
component of the CDM. Furthermore, it is found that the PBHs with quite large masses 
(~ 1O 5 M ) and very narrow mass spectrum can be formed and these can be the seeds of 
SMBHs. 

In this paper we have derived the scalar potential in the frame work of supergravity. 
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However, we can build the model without supersymmetry if we start with the potential (3). 
The Hubble induced mass terms (5) which are necessary for generating the blue-tilted 
power spectrum can be obtained through couplings with the inflaton field. For example, 
suppose that a scalar ip causes chaotic inflation and its potential is given by V(ip) = \(p 4 . 
Then the term like gtp 2 \Q\ 2 (g: small coupling) lead to the Hubble induced mass term for 
$ if we take appropriate g. 

Acknowledgment 

We thank Fuminobu Takahashi for useful discussions. This work is supported by Grant- 
in- Aid for Scientific research from the Ministry of Education, Science, Sports, and Culture 
(MEXT), Japan, No. 14102004 (M.K.), No. 21111006 (M.K.) and also by World Premier 
International Research Center Initiative (WPI Initiative), MEXT, Japan. N.K. is sup- 
ported by the Japan Society for the Promotion of Science (JSPS). 

References 

[1] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
[arXiv:1001.4538 [astro-ph.CO]]. 

[2] For a review, see S. P. Martin, arXiv:hep-ph/9709356. 

[3] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[4] Ya. B. Zel'dovich and I. D. Novikov, Sov. Astron. A. J. 10 (1967) 602, B. J. Carr 
and S. W. Hawking, Mon. Not. Roy. Astron. Soc. 168, 399 (1974). 

[5] B. J. Carr, Astrophys. J. 201, 1 (1975). 

[6] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [Erratum-ibid. 46, 206 (1976)]. 

[7] B. J. Carr, K. Kohri, Y. Sendouda and J. 'i. Yokoyama, Phys. Rev. D 81, 104019 
(2010) [arXiv:0912.5297 [astro-ph.CO]]. 

[8] M. Kawasaki, N. Sugiyama and T. Yanagida, Phys. Rev. D 57, 6050 (1998) [hep- 
ph/9710259], M. Kawasaki and T. Yanagida, Phys. Rev. D 59, 043512 (1999) [hep- 
ph/9807544]. 

[9] J. 'i. Yokoyama, Phys. Rev. D 58, 083510 (1998) [astro-ph/9802357]. 

[10] T. Kawaguchi, M. Kawasaki, T. Takayama, M. Yamaguchi and J. 'i. Yokoyama, Mon. 
Not. Roy. Astron. Soc. 388, 1426 (2008) [arXiv:0711.3886 [astro-ph]]. 

[11] P. H. Frampton, M. Kawasaki, F. Takahashi and T. T. Yanagida, JCAP 1004, 023 
(2010) [arXiv: 1001. 2308 [hep-ph]]. 



15 



[12] K. Kohri, D. H. Lyth and A. Melchiorri, JCAP 0804, 038 (2008) [arXiv:071 1.5006 
[hep-ph]]. 

[13] M. Drees and E. Erfani, JCAP 1104, 005 (2011) [arXiv: 1102.2340 [hep-ph]], M. Drees 
and E. Erfani, JCAP 1201, 035 (2012) [arXiv: 11 10.6052 [astro-ph.CO]]. 

[14] J. Kormendy and D. Richstone, Ann. Rev. Astron. Astrophys. 33, 581 (1995). 

[15] J. Magorrian, S. Tremaine, D. Richstone, R. Bender, G. Bower, A. Dressier, 
S. M. Faber and K. Gebhardt et al, Astron. J. 115, 2285 (1998) [astro-ph/9708072]. 

[16] C. J. Willott, R. J. McLure and M. J. Jarvis, Astrophys. J. 587, L15 (2003) [astro- 
ph/0303062]. 

[17] N. Duechting, Phys. Rev. D 70, 064015 (2004) [astro-ph/0406260]. 

[18] M. Kawasaki, A. Kusenko and T. T. Yanagida, Phys. Lett. B 711, 1 (2012) 
[arXiv: 1202.3848 [astro-ph.CO]]. 

[19] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) [hep-ph/0 110002], D. H. Lyth, 
C. Ungarelli and D. Wands, Phys. Rev. D 67, 023503 (2003) [astro-ph/0208055]. 

[20] J. Yokoyama, Astron. Astrophys. 318, 673 (1997) [astro-ph/9509027]. 

[21] S. Kasuya and M. Kawasaki, Phys. Rev. D 80, 023516 (2009) [arXiv:0904.3800 [astro- 
ph.CO]]. 

[22] S. Kasuya, M. Kawasaki and T. Yanagida, Phys. Lett. B 409, 94 (1997) [hep- 
ph/9608405]. 

[23] M. Dine, L. Randall and S. D. Thomas, Nucl. Phys. B 458, 291 (1996) [hep- 
ph/9507453]. 

[24] A. R. Liddle and D. H. Lyth, "Cosmological Inflation and Large-scale Structure," 
(Cambridge University Press, Cambridge, England, 2000). 

[25] J. C. Niemeyer and K. Jedamzik, Phys. Rev. Lett. 80, 5481 (1998) [astro- 
ph/9709072], J. C. Niemeyer and K. Jedamzik, Phys. Rev. D 59, 124013 (1999) 
[astro-ph/9901292]. 

[26] A. M. Green and A. R. Liddle, Phys. Rev. D 60, 063509 (1999) [astro-ph/9901268]. 
[27] M. Shibata and M. Sasaki, Phys. Rev. D 60, 084002 (1999) [gr-qc/9905064]. 
[28] A. M. Green and A. R. Liddle, Phys. Rev. D 56, 6166 (1997) [astro-ph/9704251]. 
[29] H. I. Kim and C. H. Lee, Phys. Rev. D 54, 6001 (1996). 

[30] R. Bean and J. Magueijo, Phys. Rev. D 66, 063505 (2002) [astro-ph/0204486]. 



16 



T. Bringmann, P. Scott and Y. Akrami, Phys. Rev. D 85, 125027 (2012) 
[arXiv: 1110.2484 [astro-ph.CO]]. 



17 



